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A remarkable periodic solution 
of the three-body problem 
in the case of equal masses 

By Alain Chenciner and Richard Montgomery 

Dedicated to Don Saari for his (censored) birthday 
Abstract 

Using a variational method, we exhibit a surprisingly simple periodic orbit 
for the newtonian problem of three equal masses in the plane. The orbit has 
zero angular momentum and a very rich symmetry pattern. Its most surprising 
feature is that the three bodies chase each other around a fixed eight-shaped 
curve. Setting aside collinear motions, the only other known motion along 
a fixed curve in the inertial plane is the "Lagrange relative equilibrium" in 
which the three bodies form a rigid equilateral triangle which rotates at con- 
stant angular velocity within its circumscribing circle. Our orbit visits in turns 
every "Euler configuration" in which one of the bodies sits at the midpoint 
of the segment defined by the other two (Figure 1). Numerical computations 




Figure 1 (Initial conditions computed by Carles Simo) 

a:i=-a:2=0.97000436-0.24308753j,a:3=0; y=i3=-2ii=-2a:2=-0.93240737-0.86473146i 
T=12T=6.32591398, i'(0)=2, mi=m2=m3=l 
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by Carles Simo, to be published elsewhere, indicate that the orbit is "stable" 
(i.e. completely elliptic with torsion). Moreover, they show that the moment 
of inertia I{t) with respect to the center of mass and the potential U{t) as 
functions of time are almost constant. 

1. The setting 

We consider three bodies of unit mass in the Euclidean plane = C, 
with Newtonian attraction. After using Galilean invariance to fix the center of 
mass, the phase space becomes the tangent bundle to the configuration space 
X = X\{ collisions} , 

= |x = (xi,x2,x3) G (c)^ Yl^i = o| • 

We endow X with the Hermitian mass scalar product. In the case of equal unit 
masses this is 

3 

{x, y) = ^Xj - yj = x-y + iuj{x, y). 

Its real and imaginary part are respectively the mass scalar product and the 
mass symplectic structure. The isometry group 0(2) of = C acts diagonally 
on X: the symmetry S with respect to the first coordinate axis and the rotation 
Re of angle act respectively as 

S ■ (xi,X2,X3) = (xi,X2,X3), 

Re ■ (xi, X2, xa) = (e^^xi, e^^X2, e'^x^). 

The phase space will be identified with the Cartesian product X x X with 
elements written {x,y). We define the following 0(2)-invariant functions on 
phase space: 

J = x-x, J = x-y, K = yy, U = U{x), H = - U, L = + U. 

These are the moment of inertia with respect to the center of mass, half its deriva- 
tive with respect to time, twice the kinetic energy in a frame attached to the 
center of mass, the potential function, the total energy, and the Lagrangian. 
The size r = of the configurations is the norm on X. The "force function" 
U, the negative of the potential energy, is defined by 

U = 1 1 , where ri-j = \xi — Xj\. 

ri2 ri3 r23 

(We choose units so that the gravitational constant is 1.) 
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The central configurations play a basic role in this paper. They are the 
only configurations which admit homothetic motions, that is motions for which 
the configuration collapses homothetically to its center of mass. They admit 
more generally homographic motions where each body has a similar Keplerian 
motion and in particular relative equilibrium motions where the bodies rotate 
rigidly and uniformly around their center of mass. These configurations are the 
critical points of the scaled potential function U = ^/lU . Upon normalization 
these are the critical points of the restriction f7|/=i of the potential function 
to the sphere / = 1. In the case of three bodies, the central configurations 
are completely known thanks to the works of Euler and Lagrange. There are 
three collinear configurations Ei, E2, E^, distinguished by which mass sits in 
between the two others (it sits at the midpoint in the case of equal masses), 
and two equilateral ones L~^,L~ distinguished only by their orientation. 

In naming the central configurations we have already formed the quotient 
of the configuration space by the rotation group SO (2). It is well-known ([6], 
[7]) that after reduction by direct isometrics (translations and rotations) the 
three-body problem in the plane has a configuration space which is homeomor- 
phic to . This reduced configuration space - the space of oriented triangles 
in the plane up to translation and rotation - is endowed with a metric induced 
from the mass metric on configuration space which makes it a cone over a 
round 2-sphere of radius ^. Points of this sphere, henceforth called the shape 
sphere, represent oriented triangles whose moment of inertia / is 1. This sphere 
is to be thought of as the space of oriented similarity classes of triangles [7]. 
The shape sphere is depicted in Figure 2 in the case of equal masses; it will be 




Figure 2 
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studied from a Riemannian point of view in paragraph 5. Its main features are 
the North and South poles, corresponding to the the Lagrange points and 
L~, the equator corresponding to coUincar triangles and the three meridians 
corresponding to the three types of isosceles triangles. Each meridian Mj 
intersects the equator at an Euler point Ei and an antipodal collision point Cj. 
The figure also shows the equipotential curve (level curve of U\i=i) containing 
the Euler points and a level surface of U corresponding to a higher value. 
The projection of the "eight" orbit onto the shape sphere closely resembles 
the equipotential curve passing through the three Euler points, and shares 
its symmetries. We have stressed the half-line and the meridian plane 
Ml because they explicitly appear in the proof. Note that we used a false 
perspective for the poles to increase readability. 



2. The orbit 



Let T be any positive real number. We define actions of the Klein group 
Z/2Z X Z/2Z on M/TZ and on as follows: if a and r are generators, 

T T 

a-t = t + —, T-t = -t + —, a-{x,y) = (-x, y), r • (x, y) = (x, -y). 

Theorem. There exists an ''eight" -shaped planar loop q : {R/TZ),Q 
B?,0 with the following properties: 

(i) for each t, 

q{t) + q{t + T/3) + q{t + 2T/3) = 0; 

(ii) q is equivariant with respect to the actions o/Z/2Z x Z/2Z on R/TZ and 
M? defined above: 

q{u ■ t) = a ■ q{t) and qij ■ t) = t ■ q{t); 

(iii) the loop x : R/TZ — > X defined by 

x{t) = {q{t + 2T/3), qit + T/S), q{t)) 

is a zero angular momentum T-periodic solution of the planar three-body 
problem with equal masses. 

The rest of the paper is devoted to a proof of this theorem. 



A REMARKABLE PERIODIC SOLUTION 



885 



3. The structure of the proof 



The main ingredients of the proof axe the following: 

(i) The ^^direct" method. One twelfth of the orbit is obtained by minimiz- 
ing the action (we note T = T/12) 



over the subspace A of H^{[0,T], X) consisting of paths which start at an 
Eulcr configuration, say (3 in the middle) with arbitrary size and end at an 
isosceles configuration, say Mi (ri2 = ria), again of arbitrary size. Existence 
of a minimizing path is standard. The main point is to show that such a path 
has no collision. 

(ii) Reduction. The kinetic energy can be expressed as the sum of two 
nonnegative terms: K = ifj-ed + ^rot- (This is Saari's decomposition of the 
velocity. See paragraph 5 below.) ii'rcd corresponds to a Riemannian metric 
on the quotient space X /S0{2) induced from the metric K on X. It is the 
deformation part of the kinetic energy, including the homothetic part. (See 
Lemma 5 below for its explicit expression.) K^ot is the rotational part of 
kinetic energy. Because we consider the planar problem, Kj-ot = with 
C the angular momentum vector. (This equality is replaced by the inequality 
Krot > I^P/-^ ill the spatial three-body problem. See Sundman's inequality 
IK — > in [1].) The boundary conditions defining A are invariant under 
rotation, which means that if x{t) is such a path, so is g{t)x{t) for g{t) any 
sufficiently regular curve in S0(2). Changing the choice of g{t) changes i^rot 
while leaving Kred fixed. By appropriate choice of g{t) we can insure that 
^rot = and hence C = 0. This shows that any minimizer for the problem 
of (i) has zero angular momentum. Moreover, the original problem has been 
reduced to the problem of minimizing the reduced action 



over paths lying in the reduced configuration space A'/S0(2), and satisfying the 
same boundary conditions. The advantage of this reduction is that the reduced 
configuration space is topologically M^, hence paths are easy to visualize (see 
Figures 2, 3, and 4). 

(iii) Comparison with Kepler to exclude collisions. Instead of computing 
local variations of the action, we compute the infimum of the action over all 
paths in ^ ( [0, T] , A") with collision. Denote this infimum by A2 , the "2" stand- 
ing for two bodies because we explicitly compute A2 via a two body problem. 
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We compare A2 to the action a of an explicit carefully chosen collision-free 
"test path" in A. We show that A2 > a, yielding the result that the minimizer 
must be collision-free. 

(iv) Symmetries and area rule. The equality of masses gives a supply 
of symmetries corresponding to interchanging the masses. Using these, we 
construct eleven other congruent copies of the minimizer from (i). The first 
variation of action formula shows that these copies fit together smoothly with 
the original to form a single orbit, periodic in the reduced configuration space 
(i.e. modulo rotations), of period 12T. To reconstruct the motion in the inertial 
plane, i.e. in X, we use the symmetries combined with the "area formula" (see 
[9]). These tools yield that the motion in X is periodic, that all three masses 
indeed move along one and the same curve in the inertial plane, with the Klein 
group symmetry described in the theorem. 

(v) Proving the curve has the shape of a figure eight". We prove that 
the individual angular momentum of any one of the three bodies is zero only 
as this body passes through the origin. This implies that each one of the two 
lobes of the curve is starshaped. 

Other proofs. We have presented the shortest proof we know. Another 
proof proceeds by local perturbations so as to destroy collisions while decreas- 
ing the action, getting rid of the binary and triple collisions separately. Special 
care must be taken with the direction of perturbation used when destroying 
binary collisions. Estimates arc significantly longer than those given here, but 
do not require numerical integration. On the other hand, they do not yield the 
estimates on the orbit given in Appendix 1. 



Let us recall that A is the subspace of H^{[0,T],X) consisting of paths 
which start at the Euler configuration E3 with arbitrary size and end at an 
isosceles configuration of type Mi with arbitrary size. This paragraph and the 
next one are devoted to proving 

Proposition. A path in A which minimizes the action 



4. The exclusion of collisions 




has no collisions. 



Proof of the proposition. Surprisingly, we are able to treat double and 
triple collisions simultaneously. The first key point is the (trivial) remark that 



A REMARKABLE PERIODIC SOLUTION 



887 



the action 



A{mi,m2,m3;x) = J {^^rni\xi{t)\'^ + ^ 



rriimj 



,Xi(t) — Xj(t)\ 
l<i<j<3 ' ■' ''I 



dt 



along a given path 1 1— > x{t) = [xi{t),X2{t),X3{t)) is an increasing function of 
any of the masses. Setting, for example, mi to zero, yields 

A{x) = Ail, 1, 1; a;) > ^(0, 1, 1; x). 

The last term is the action of a 2-body problem with equal masses. We use 
this remark in the following lemma: 

Lemma 1. If x E H^{[0,T], X) has a collision, either double or triple, 
then its action is greater than A2, where A2 is the action of the two-body 
coUinear solution in which two unit masses start at collision and end with zero 
velocity at a time T later, the center of mass being fixed throughout. ( This 
might be called half a collision-ejection elliptic orbit.) See below for an explicit 
formula for A2. 

Remark. The lemma asserts that the infimum of the action A{x) over 
all collision paths x G H^{[0,T],X) is greater than or equal to A2. Indeed 
this infimum equals A2- Imagine a sequence x„ of paths in which m.2 and 
m3 perform half the Keplerian collision-ejection orbit described in the lemma, 
while mi remains fixed, and at distance n from the 2-3 center of mass. We 
have A{xn) A2 as n ^ 00. 

Proof of Lemma 1. Let us suppose that masses 2 and 3 (and possibly 1) 
collide at the instant Ti. As just explained, we lower the path's action by 
setting mi equal to zero. We can forget about the position of mass 1 for this 
new action, and are left with the action for the Keplerian two body problem, 
investigated by Gordon in [4]. According to Gordon each part of the curve x, 
before Ti and after Ti, has action greater than or equal to the corresponding 
collinear motion of masses m2 and in which they collide at Ti and are at 
rest at their other endpoints, t = or t = T. Indeed, by doubling each part by 
concatenating it with itself but reversed, we get two closed paths, each going 
from collision to collision, one in time 2Ti and the other in time 2(T — Ti). The 
absolute minimum for the collision-to-collision problem was shown by Gordon 
to be the collision-ejection solution. Its action is proportional to T^/^, a convex 
function of the period T. Following Gordon, this convexity implies that the 
action is lowered further if we replace our previous two motions by the single 
collision-ejection solution beginning and ending at collision with no collisions 
in between. Half of this path realizes the infimum to collision in time T for 
the Kepler problem. This ends the proof. 
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The next lemma introduces our cohisionless test path and reduces the 
proof of the proposition to an estimation of the length Iq of its projection on 
the shape sphere. 

Equipotential test paths. Fix I = Iq and U = Uq where Uq is the value 
of the potential function at any one of the Euler configurations on the sphere 
I = Iq. Viewed in the reduced configuration space, this defines a curve on the 
two-sphere of radius y/lQ (see Figure 2). Take the one-twelfth of this curve 
lying above the equator, connecting £'3 to Mi. Traverse this curve at constant 
speed, the speed chosen so as to finish at the desired time T. This gives us a 
family of reduced test paths in A'/S0(2) depending on Iq. The corresponding 
paths in X are those which have zero-angular momentum and project to these. 
The lengths of these paths are £oV^ where Iq is the length of the path when 
/o = 1, henceforth called the "Euler equipotential length." 

Lemma 2. The minimum a of the actions of equipotential test paths is 
smaller than A2, the infimum of the actions of collision paths in time T, if and 
only if the Euler equipotential length Iq satisfies 

eo<l 
5 

Proof of Lemma 2. The action A2 of Lemma 1 is half that of the collision- 
ejection path of period 2T. It is computed (for example in [2] ; take k = 
there) to be: 

A2 = ^x3x(27r2)5(lc/2)i(2r)l; = ^ . 

The constant U2, which is called Uq in [2], is the (constant) value of the scaled 
potential function U = ^/lU for the 2-body problem with both masses equal 
to L 

We now evaluate the minimum a of the actions of the equipotential test 
paths. (This computation is the same as the one in [2].) In computing the 
action A{Iq) of the test path at radius \/Iq note that both integrands are 
constant. The action is then 

A{Iq) = {^Ko + Uq)T, where Kq = (^-^) , Uq = 



The constant = is the value of U = y/lU at the Euler configurations. 
As in [2], we are left with the problem 

... .... 1 ( Ig^/Tq \\ , 5 1 

mmimize:^(/o) = -(^^J r + -— T. 
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This function has a unique minimum with respect to /q at 



2 

3 

4 



The corresponding action is 
Finally, a < ^2 if and only if 
which amounts to 



5. Length computations 

Lemma 3. The Euler equipotential length satisfies io < ir/5. 

In order to get this estimate we will need explicit coordinates on the 
quotient space X/S0{2), and expressions for the metric and the potential 
function in these coordinates. The estimation of Iq is then obtained numerically 
with great accuracy. 

(i) The quotient map. One realizes the quotient by composing the "Ja- 
cob! map" J with the "Hopf map" /C. The configuration space X is a two- 
dimensional complex Hermitian vector space. As such, it is isometric to C^. 
Jacobi coordinates 

J -.X 

defined by 

{Z1,Z2) = J{X1,X2,X3) = ^"^N - X2), (^Xl - 

realize this isomorphism. (Recall that if (xi,X2,X3) is a point in X then xi + 
X2+ X3 = 0.) Being an isometry, we have / = |2;ip+|2:2p in Jacobi coordinates. 
The action of SO (2) corresponds to the diagonal action of the complex unit 
scalars on Jacobi coordinates {zi,Z2). That is to say, X/S0{2) = C'^/S^. 
The quotient by rotations is realized by making a vector out of the invariant 
polynomials for this action: 

}C{zi,Z2) = {ui,U2 + iu3) = - 1 2:2 2^12:2). 
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This is the Hopf (also called Kustaanhcimo-Stiefel) map 

Remark. A compact definition of /C is obtained by identifying (re- 
spectively R^) with the quaternions H (respectively with the purely imaginary 
quaternions H) as follows: 

3 {zi,Z2) ^->■ q = zi + Z2j E B. 
3 {ui,U2,us) I— s- mi + {u2 + iu^)k = ui — u^j + U2k € H. 

Then /C : H ^ H is defined by /C(g) = qiq. 

Here are some properties of these mappings: 

|/C(Z1, Z2)|' = ul + ul+ul = (|Z1|2 + \z2\''f = I^. 

The 3-sphere I = 1 is sent by /Co J" to the unit 2-sphere of (the shape sphere) 
according to the Hopf fibration. Indeed, restricted to / = 1, our formula for /C 
is the standard formula for the Hopf fibration. The location on this sphere of 
the collision points, the Euler points and the Lagrange points are: 

Ci = (-1,0,0), C2 = (^,^,0), C3 = (^,-^,0), 

Si = (1,0,0), = (-1,-^,0), = (-1,^,0), 
L+ = (0,0,1), L- = (0,0,-1). 
Using the above formulas and the expressions 

r23 = V2\zi\, rai = \^/3/2z2 + (1/^2)^1 1, = \VV2z2 - {l/V2)zi\, 

the proof of the following lemma is immediate: 

Lemma 4 (Hsiang [5]). To u = {ui,U2,U2) in the shape sphere, cor- 
responds a triangle with sides r23 = \/l — Ci • u, r^i = ^1 — C2 ■ u, ri2 = 
^/l — Cs ■ u, where the scalar product is the standard Euclidean one in R^. 

(ii) The orbit metric. Wc derive a formula for the reduced metric i^'icd by 
computing the distance d{x,y) between the SO(2)-orbits of x and y in X. As 
S0(2) acts by isometrics, 

d^{x,y) = mi "^{xi-e'^yif 

i 

= inf + |yp — 2x ■ j/cos 9 + 2lv{x, y) sin 9\ . 
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Taking the derivative with respect to 0, we see that the minimum occurs 
sA, 6 = 9o where a; • y sin + ^{x, y) cos 6q = 0. This imphes that 



(f{x,y) = + - 2^/{x ■ yf+u:{x,yf = \x\'^ + \y\'^ - 2|(x,y)|. 

The term in the expansion of d?{x,x + sv) is the reduced kinetic energy 
Kj.gd{x, v) corresponding to the decomposition i^j-ed = K — Kj-ot- We find 

(This is the expression of the pullback of the natural induced metric on the 
quotient X/S0{2) by the projection of X onto the quotient.) This is consis- 
tent since Kj-ed{x,v) = when v is tangent to the orbit of x, i.e. when v is 
proportional to ix. 

(iii) The length Iq in spherical coordinates. It will be convenient to use 
spherical coordinates defined in the shape space R'^ by 

ui = cos (p cos 6, U2 = cos sin 0, = sin ip. 

We have 

which justifies the choice of the notation r. A tedious but simple calculation, 
or an appeal to the ?7(2)-invariance of i^red) now proves 

Lemma 5 (see [7]). In spherical coordinates the quotient metric corre- 
sponding to the reduced kinetic energy -R'red occurring in the reduced action is 
given by 

2 

ds"^ = dr"^ + ^{cos^ipde'^ + (i<y?^). 

In particular the shape sphere I = r^ = \ is isometric to the standard sphere 
of radius 1/2, and the shape space is the cone over this sphere., the sphere 
itself consisting of all those points at distance 1 from triple collision. 

Using Lemma 4 we can express the equipotential curve through the Euler 
configurations on the shape sphere by the implicit equation: 

1 1 

+ 



VI + cos cos 6/ ^i + cosv^cos(e + f ) 



l + cos^cos(^ + ^) v2 

This curve is a double covering of the equator 99 = and as such may be 
parametrized by a function = (^(6), provided is allowed to vary in an 
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interval of length 47r. We are interested in io which is one twelfth of its length. 
It follows from Lemma 5 that 

2i7V — 

4 = ^/ Vcos2 ip{e) + ^^e) de = l- [ ' ^ycos^ ^{9) + ^^e) de. 

1^ Jo ^ Jo 

To finish the proof of Lemma 3, and hence of the existence of the collision 
free minimizer, we use the following numerical estimate of the Euler equipoten- 
tial length £o obtained by Carles Simo and later confirmed by Jacques Laskar: 

<io< . 

5.082553924511 - " - 5.082553924509 

These estimates were obtained by using a Newton method for computing (p{d) 
and ip'{0), and then the trapezoid method for computing the integral. 

Remark. We explain the meaning of the spherical coordinates in terms of 
triangles. The parallels or "latitudes" (p = constant in the shape sphere cor- 
respond to triangles with the same orientation and the same ellipse of inertia 
up to rotation. Indeed, this set of triangles is characterized by a common area 
(see [1]). But the area is proportional to lmziZ2, that is to U3 = sini^, the 
height function on the sphere. The meridians or "longitudes" 6 = constant in 
the shape sphere are defined by a linear relation between the squares of the 
mutual distances. These properties of the coordinates {0, tp) are a consequence 
of the invariance of the metric under the orthogonal group 0(P) of the dispo- 
sition space T) (see [1] for the definition). The equilateral triangles L± (north 
and south poles) are fixed points of the action of SO(P) and the parallels above 
are the circles with center L± which are also the orbits of the action of SO(P). 
The geodesies through L± orthogonal to these circles are the meridians. They 
are transitively transformed into each other by SO(X') and each one is fixed by 
an involution in 0(P). 

6. Symmetries: Proof of the existence of the "eight" 

Recah our action of the Klein group Z/2Z x Z/2Z on M/TZ and M?. li a 
and r are generators, 

T T 

a - t = t + T -t = -t+ —, a-{x,y) = {-x,y), r ■ {x,y) = {x,-y). 

Lemma 6. After being symmetrized according to the pattern of the Euler 

equipotential curve on the shape sphere, a minimizing path x gives a T = 12T 
periodic loop {still called x), with zero angular momentum, which, up to a time 
translation and a space rotation, is of the form 

x{t) = {q{t + 2r/3), q{t + T/3),q{t)) , 

described in the theorem. 
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The proof of Lemma 6 proceeds in three steps. 

Step 1. Observe that the minimizing arc is orthogonal to the two man- 
ifolds Es and Mi constraining its endpoints. This follows from the boundary 
term arising in the first variation formula for the action. 

Step 2. Observe that upon reflecting this arc about one of the three 
meridians, or about the equator, we will obtain another minimizing solution 
arc, one with permuted endpoint conditions; e.g. with Ej and in place of 
E3 and Ml. Using these reflections we build the entire closed solution curve in 
the reduced configuration space. It consists of 12 subarcs all congruent to our 
original minimizcr. Orthogonality guarantees that they fit together smoothly, 
thus forming a single solution. 

More precisely, because reflection about the meridian Mi is a symmetry of 
the reduced action (and hence of the equations), and because the minimizing 
arc is orthogonal to the meridian at its endpoint, when we continue the solution 
represented by the arc through the meridian Mi , the result is the same as if we 
had reflected it about the meridian, and then reversed the direction of time. 
In symbols: x{T/12 + t) = s\{x{T/Vl — t)) where si is reflection about the 
meridian Mi, and where T = 12T will be the period of the full orbit. (T is the 
time it takes to hit the meridian.) 

The reflection si can be realized in the inertial plane as follows: at time 
T = T/12 the triangle is an isosceles triangle of type Mi and hence has a 
reflectional symmetry r. Choose coordinates in so that t(x, y) = {x,—y), 
i.e. so that the perpendicular bisector of the edge joining 2 to 3 is the x-axis. 
Let 5'i(xi, ,T2, .X3) = {xi,X3, X2) be the operation of interchanging masses 2 and 
3. Then si = Si o t = t o Si. 

We now have a solution from E^ to E2 in time 2T = T/6. By a similar 
argument, to continue this arc of solution through E2 we must perform a 
half-twist H2 through E2 and reverse time: x{2T - t) = H2{x{2T + t)). This 
half-twist is a symmetry of the action being the composition of reflection about 
the equator with reflection about the meridian M2. It is realized in the inertial 
plane by interchanging masses 1 and 3 and then performing the inertial half 
twist a o t(x, y) = (— x, —y) about the origin. 

We obtain in this way an arc of solution from E^ to Ei in time 4T = T/3. 

Continuing around the equator in this manner with the appropriate reflec- 
tions or half twists we construct a smooth curve in the reduced configuration 
space which consists of 12 congruent arcs, alternating in pairs above and below 
the equator, so as to have the same symmetry as the equipotential curve. It is 
a solution curve, and is T-periodic mod rotations. 

Step 3. We have constructed the projection of our solution curve to the 
reduced configuration space. We now reconstruct the full solution curve, show 
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that it is periodic (i.e. in inertial space), and show that it satisfies ah the 
properties described in the theorem. This is done by invoking the area rule for 
reconstructing the original dynamics from the reduced dynamics, and by using 
the symmetries of the curve. 

Figure 3 shows segments of the reduced orbit on the shape sphere and 
anticipates the reconstructed orbit in the inertial plane. 




3 



Figure 3 

The area rule tells us how to recover the motion of the masses in the 
inertial plane, given the curve representing this motion in the shape sphere. 
Suppose the shape curve is closed. Then the initial and final triangles in the 
plane are similar. The angle of rotation which relates these two triangles up 
to scale equals twice the spherical area enclosed by the shape curve. (The area 
of the sphere of radius 1/2 is vr. The factor of 2 in the area formula insures 
that the answer is well-defined modulo 27r.) For a proof of the area rule see, 
for example, [9] or references therein. If instead the shape curve is not closed, 
but begins and ends on the equator of collinear shapes, then we close it up by 
travelling "backwards" along the equator from the endpoint until we reach the 
beginning point. Now compute twice the signed area enclosed by this closed 
curve. This equals the angle between the two lines in the inertial plane which 
contain the initial and final configuration for any zero angular momentum curve 
realizing the given shape curve*. Finally, if the curve begins or ends on one 
of our three isosceles meridians, we compute the angle between the beginning 



* There are two ways to close up the curve into a loop, depending on which way we travel the 
equator. The two angles so computed differ by vr, this being twice the area of a hemisphere of a 
sphere of radius one-half. This is no problem since the angle between two unoriented lines is only 
defined modulo tt. 
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and ending symmetry axis of the isosceles triangle by following the appropriate 
meridian up or down to the equator, travelling along the equator to close up 
the curve, and then computing the area within the resulting closed curve. 

The fact that the signed areas depicted on this figure are equal to zero 
implies that as we travel our solution curve 

(i) if we start at an Euler configuration and follow the orbit for a time T/3 = 
4T, passing through an intermediate Euler configuration at time 2T we 
arrive at an Euler configuration with the three masses sitting on the same 
line as that of the initial Euler configuration. That is to say, there is no 
rotation of the Euler line, contrary to what happened at the intermediate 
time 2T. 

(ii) after time r/2, an isosceles configuration returns to itself, up to reflection 
(equivalently up to interchange of the symmetric vertices). There is no 
rotation of the symmetry axis of the triangle. 

Choose the origin of time t = to correspond to being in the Euler 
configuration E-^. Set q{t) = x-i,{t) where x{t) = {xi{t),x-2{t),xz{t)) is our 
solution. The first property implies that 

q{t) = xs{t) for < t < T/3, 
q{t) = X2{t - T/3) for T/3<t< 2T/3, 
q{t) =xi{t- 2T/3) for 2T/3 <t<T. 

Thus, after time T/3 (respectively, 2T/3), the bodies 2,3,1 have been 
replaced by the bodies 3,1,2 (respectively, 1,2,3) with the same velocities. 
The three bodies move along the same closed curve q{t) of period T with a 
phase shift relative to each other of T/3. Using (ii), the Klein symmetry follows 
easily. 

Figure 4 shows the projection (still called x) of the orbit in the reduced 
configuration space. 

Step 4. It remains to be proven that the equivariant curve q we con- 
structed not only has the required symmetry but also has the shape of a figure 
eight without extra small loops or other unpleasant features. We will use the 
following basic lemma: 

Lemma 7. The angular momentum q{t) A q{t) is nowhere vanishing for 
t < < T/4; i.e. in any one- quarter of the curve the angular momentum of 
the mass tracing out that quarter is nonzero. 
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x{-t + 27/6) 
x{-t + 2T/3) 




t + 5T/6) 



5T/6) 



Figure 4 

Proof of Lemma 7. We will need the fact that Newton's equations are 
satisfied, and two consequences of the minimality of the curve in the shape 
sphere. 

We begin by computing the time derivative of g A g. It is 

d (Pq 

Now use the fact that q{t) = xs{t) where {xi{t),X2{t),X3(t)) satisfy Newton's 
equations. Also use the fact that the center of mass is zero at all times: 
xi + X2 + xs = 0. These yield 



It 



1 



1 



(X3 A Xi). 



13 ' 23 - 

There are only two ways this quantity can be zero: either 

(A) xi and X3 are linearly dependent , or 

(B) ri3 = r23. 

To eliminate these possibilities, we use the refiection principle in shape space. 
First, if xi and X3 are linearly independent, then the entire configuration is 
collinear. Thus it crosses the equator in shape space. This can happen for 
a minimizing arc only at an Euler point. For if it did at another point, the 
minimizing arc between Euler and Isosceles (0 < t < T/12) would be divided 
into two (or more) subarcs, one of which lies below the equator and the other 
above. Refiect one of these arcs across the equator, leaving the other(s) fixed. 
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The resulting arc has the same action, has no collisions, and still connects Euler 
to Isosceles, and hence is also a collision-free minimizer. But it is no longer 
analytic, contradicting the fact that collision-free minimizers correspond to 
solutions. 

Next suppose that ri3{t) = r23{t) at some time t, with < t < Tj 12. 
This asserts that the curve has recrossed the isosceles meridian passing through 
our initial Euler point. The same reflection principle applies, with the Euler 
meridian playing the role of the equator. 

Wc can now suppose that, as in Figure 4, our minimizing arc lies in the 
left upper quarter of the shape sphere for < t < r/6. The arguments thus 
far show that the angular momentum q{t) A q{t) decreases from the value 
for < t < T/6 and increases for T/6 < t < T/4 (a bivector in the plane is 
called positive if it is a positive multiple of the standard area form ei A 62). It 
remains to notice that it takes a negative value at i = T/4 to conclude that it 
stays strictly negative for < t < T/4. It follows that it stays strictly negative 
for < t < T/2 and strictly positive for T/2 <t<T. 

Corollary. Each lobe of the curve is starshaped: the only time Xi A Xi 
becomes zero, for i = 1,2, 3, is when xi passes through the origin. 

Proof of the corollary. In polar coordinates (r, 6) the angular momentum 
has the well-known expression x Ax = (r^^)ei A 62. From this it follows that 
the polar angle 9{t) of the curve q{t) decreases monotonically over the interval 
(0,T/2) from its maximum value of ^(0) when r(0) = to its minimal value 



We get estimates on U and K along our solution. These imply estimates 
on I and J by using the fact that the solution has zero mean, which follows 
from its symmetries. 



of e(T/2) 



^(0). 



Appendix 1. Estimates 



Let 



2 




be defined as above and 




a = i^Ko + Uo)T = ^UoT = -ZH^T, 



Let (/) = i f{t) dt be the mean value of a function [0, T] — > R. 
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Lemma 8. A minimizing path satisfies the following estimates: 
{U) = {K)<Uo = Ko, H>Ho, {l)<^Io, (|J|) < ^Jq. 

Proof of Lemma 8. Because the conservation of energy is true almost 
everywhere on a minimizing path*, such a path x has action 

A{x) = HT + 2{U)T <a = HqT + 2UqT. 

But we deduce from the Lagrange- Jacobi identity J = \l = 2H + U = K — U 
that 

{K) - {U) = J{T) - J(0). 

Because of the symmetry of the minimizing path, J(T) = J(0) = (note that 
the eventual presence of a double collision at the end of the path would not 
change this fact because then J{t) would behave as {T — t)^ which goes to 
as t goes to T). This implies 

{K) = {U) = -2H, 

so that 

A{x) = -3HT <a = -3HoT , hence H > Hq. 

The inequalities concerning ([/) and (iC) follow immediately. 

To bound (/), note that by construction our x has zero average (in X) 
over its full period 12T. By the Poincare lemma this implies {K) > 

and consequently we get the estimate on ^/). 

Finally, Sundman's inequality IK — > yields the bound on (| J|). 

Elmo's numerical computation of actions. The following numerical esti- 
mates of various important actions were obtained by C. Simo. The period 
here is taken to be T = 2it/\2. 

^2 = (^)^/^| = 2.0583255..., 
a = A{test) = (2257rZ§/32)V3 = 2.0359863 . . . , 
^min = ^(solution) = 2.0309938 .... 

It is interesting to notice that the action of any element of A undergoing a 
triple collision is much higher. Indeed, Sundman's inequality implies it is 
higher than the action ^3 of an equilateral homothetic solution from collision 



* This classical result can be proved by computing the variation of action due to a change of 
parametrization. 
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to zero velocity. This last action is given by the same formula as A2 with the 
value U3 = 3 oi U at the equilateral configuration replacing U2 = l/\/2, so 
that 

2 

^3 = (^V2) ' X ^2 = 5.39433 .... 



Appendix 2. How this orbit was discovered 



One of us (R.M.) had been searching for several years for periodic orbits 
in the three-body problem using the method of minimization of the action over 
well-chosen homotopy classes of loops in the configuration space. (See [7] where 
the collision problem is avoided by a strong force hypothesis, and where poten- 
tially interesting homotopy classes are described.) Approximately at the same 
time, but independently, both authors realized that equality among the masses 
could make the variational approach more tractable by allowing us to impose 
additional symmetries on competing loops. This led to the preprint [3] by A.C. 
and A. Venturelli, submitted to Celestial Mechanics, in which new periodic or- 
bits were found for the spatial four body problem with equal masses, and to the 
preprint [8], submitted by R.M. to Nonlinearity. A.C. was asked to act as a ref- 
eree for preprint [8] , titled "Figure eights with three bodies" , which described 
two different types of periodic orbits of the three-body problem according to 
whether or not the masses were all equal or only two were equal. Only the orbit 
for all equal masses survived careful scrutiny by A.C. and A. Venturelli. The 
other orbit - which, curiously, was the one which had given the paper its title 
- was supposed to be a figure eight not in the plane, but in the shape sphere. 
However the proof of the absence of collisions for this orbit was found to be 
in error. In trying to understand the case of equal masses, A.C. discovered, 
at first experimentally and then mathematically, that the three equal masses 
must travel along a fixed eight-shaped curve in the plane. This discovery gave 
new life to the title of the preprint. The numerical experiment grew out of an 
example called "figure eight attractor" in the nice program "Gravitation" by 
Jeff Rommereide. The success of the experiment came from the constraints 
placed on the velocities v at any Euler point of the orbit. These constraints, 
depicted on figure 1 for as vi = u 2 = —V3/2, are due to the fact that the an- 
gular momentum is zero and that each Euler configuration is an extremum of 
the moment of inertia / along the orbit: dl{v) = 0. By making more stringent 
symmetry assumptions on the orbit than R.M had made, A.C. was then able 
to give a direct and simple proof, partly in the spirit of [2], of the absence of 
triple collisions, and to obtain estimates for / and U along the orbit. Finally, 
R.M. noticed the trick of forgetting one mass, which made the calculations for 
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triple collisions extend to double collisions and bypassed completely any local 
variational analysis. Precise numerical computations by Carles Simo, using 
the special form of velocities at an Eulcr configuration, gave accurate pictures 
of the eight and showed in particular its "stability"*. 
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